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Abstract: We present finite temperature solutions describing N c D5 branes wrapped 
on the S 2 of the resolved conifold in the presence of Nf flavor branes sources and 
their backreaction i.e Nf/N c ~ 1. In these solutions the dilaton does not blow up 
at infinity but stabilizes to a finite value. Thus, we can use them to generate new 
ones with D5 and D3 charge. The resulting backgrounds are non-extremal versions of 
the "flavored" resolved deformed conifold. It is tempting to interpret these solutions 
as gravity duals of finite temperature field theories exhibiting non-trivial phenomena 
as Seiberg dualities, Higgsing and confinement. However, a first necessary step in this 
direction is to investigate their stability. We study the specific heat of these new flavored 
backgrounds and find that they are thermodynamically unstable. Our results on the 
stability also apply to some of the non-extremal backgrounds with Klebanov-Strassler 
asymptotics found in the literature. 
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Appendix: General form of B 2 



1 Introduction 



The gauge/gravity duality [1],[2],[3] provides a completely new framework for the un- 
derstanding of strongly coupled field theories. It has opened a new area of study at 
the interface of field theory and string theory. The applications are numerous and 
include formal developments as well as phenomenological topics. One of the goals is 
to understand aspects of QCD, at zero and finite temperature, through the use of a 
suitable string dual. In this spirit, the two landmarks for non-conformal backgrounds 
with Af = 1 supersymmetry — the deformed conifold model of Klebanov and Strassler 
(KS) [4] and the wrapped branes model of Maldacena and Nunez (MN) [5] - have 
been thoroughly studied and the connections between them well understood. Both 
backgrounds, MN and and KS, are special cases of a general ansatz proposed by Pa- 
padopoulos and Tseytlin in [6]. In [7], a family of solutions interpolating between KS 
and MN was found by imposing supersymmetry conditions on the ansatz of [6] . In the 
field theory side this family of solutions was shown to correspond to different vacuum 
expectation values of a baryonic operator. Varying the string coupling constant of this 
family of SUGRA solutions smoothly interpolates between KS and MN. In order to 
incorporate dynamical flavor in the fundamental representation in these backgrounds 
we have to consider the backreaction of flavor branes. This difficult problem can be 
tackled using a smearing technique and has been extensively explored in the literature 
[8-14]. 

In [15] the authors presented a particular solution of wrapped D5 branes that inter- 
polates between the deformed conifold with flux and the resolved conifold with branes. 
Applying U-dualities to this solution adds D3 charge and we recover the Klebanov- 
Strassler baryonic branch. Thus, the chain of dualities takes us from a background 
with only dilaton and three form flux Ht s \, to one with dilaton, H^, Fr 3 ) and F/$\. 
This procedure -also refer to as rotation x - requires that we start with a D5 brane 
background whose gravity modes and field theory modes are not decoupled. Unlike 
MN and the flavored deformations of MN [8], this type of solutions have a stabilized 
dilaton. In this framework, the decoupling limit is taken after the dualities have been 
applied so it is only the backgrounds with dilaton, H@), F( 3 ) and that will have an 
interpretation as field theory duals [15]. This solution generating technique can also be 
used to study backgrounds with flavor and backreaction. In [16] the authors presented 
a new flavored wrapped brane solution with stabilized dilaton. Using this new solution 
as a seed for the rotation procedure, a flavored generalization of the KS baryonic branch 
can be obtained. The dual field theory exhibits Seiberg dualities and Higgsing and was 
conjectured to be a mesonic branch of KS. The SU(Nf/2) flavor symmetry arises only 

1 As shown in [16] it is a rotation in the space of Killing spinors. 



- 3 - 



in the IR, at the bottom of the cascade after Seiberg dualities and Higgsing have taken 
place; the quarks are really bi-fundamentals but since one of the gauge groups is very 
weakly coupled it can be thought of as flavor. This scenario was further developed in 
[17],[18],[19]. 

The rotation procedure can also be applied to non-supersymmetric [20] or non- 
extremal backgrounds [21]. Thus, it provides an alternative way of generating non- 
extremal deformations of KS which is a fascinating subject with rich new physics [22]. 
In particular, KS black holes were studied in [23], [24]. 

In this paper we present solutions that are non-extremal deformations of the fla- 
vored warped deformed conifold of [16], i.e we present new non-extremal flavored solu- 
tions that through a chain of dualities take us to a non-extremal cascading theory. 
We also study their stability. The backgrounds presented here contain the finite- 
temperature solutions of [21] as a special case (Nf = 0, a(r) = 0) 2 . On the other 
hand, in our solutions the non-extremal parameter enters in the warp factor, Fr 3 \, 
and F( 5 ) and thus, they are not (for Nf = 0) of the type studied in [23], [24]. 

We study the specific heat of these new backgrounds and find it to be negative 
both before and after the rotation. In [25] the authors showed that non-extremal MN 
is unstable and in [26] this instability was related to the existence of a tachyonic quasi- 
normal mode . Even though in this work we consider different wrapped D5 solutions 
than the ones studied in [25], [26], it is possible that the same effect is present here. 
Also in [25] the authors argued the existence of a phase transition above which the 
theory is in a high temperature thermodinamically unstable phase. It is possible that 
the U-dualities map the phase transition of [25] to a phase transition in a cascading 
non-extremal background [27], [22], [26], [23], [24]. We leave these issues for future study. 

The remainder of the paper is organized as follows. In Section 2 we review the 
general ansatz for wrapped fivebranes branes at zero temperature with and without 
flavor. This ansatz is very general and different types of solutions have been explored 
in the literature. In [8] solutions with flavor and backreaction were obtained using a 
smearing procedure. The solutions of [8] have the standard linear dilaton behavior. In 
[28] and [16] it was shown that there exist also flavored solutions with a dilaton that goes 
to a constant in the UV; these are the type of backgrounds to which the U-dualities can 
be applied and the ones that are relevant in the present work. In Section 2.2 we present 
the equations of motion for a non-extremal deformation of the flavored backgrounds of 
[16]. In Section 3 we obtain a numerical solution to the equations of Section 2.2. Once 
these solutions are under control we add D3 charge through the U-dualities and arrive 

2 a{r) is a function appearing in the metric whose infrared behavior is related to the parameter £ 
used in [7] to parametrize the baryonic branch. For small values of r, a(r) ~ 1 — £r 2 + C(r 4 ) [7]. 
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at a non-extremal flavored deformed resolved conifold. In order to numerically solve 
the equations of motion we first have to obtain a consistent UV expansion. We find a 
UV expansion with 11 free parameters that completely determine the UV behavior at 
any order. We then integrate back and match to the required regular horizon behavior. 
Having obtained the non-extremal solution we proceed in Section 4 to apply the chain of 
dualities to obtain solutions with nontrivial dialton, Fr 3 ), H^andF^. We present the 
numerical results as well as the UV behavior. In section 5 we study the thermodynamic 
stability of the seed background (wrapped fivebranes) and of the background generated 
by U-duality. We find that they are both unstable. Our solutions contain as a special 
case the ones presented in [21] thus we expect our results to hold for those backgrounds 
as well. 

2 Flavored backgrounds 

2.1 Review of flavored wrapped D5 backgrounds, T = 

In a seminal work, Casero, Nunez and Paredes [8] presented a framework to incorpo- 
rate the backreaction of the flavor D5-branes. The geometries obtained depend on the 
ratio Nf/N c which can be of order one even in the N c — > oo limit. These spaces are 
singular at the origin, but the singularity is a "good" one in the sense of the criterion 
in [29] which means that the metric component gu remains finite in the limit r — > 0. 
Everywhere else the geometry is smooth and the curvature small as long as ^A^f. 3> 1. 
These backgrounds are conjectured to be dual to Af = 1 SQCD in four dimensions 
with a large number of flavors, up to the same caveats concerning the decoupling of the 
KK modes that were already present in the discussion of the original Maldacena-Nunez 
background without flavor. 

The general strategy is the following [8]: One introduces a deformation of the Maldacena- 
Nunez background due to the presence of flavor D5-branes, derives the corresponding 
BPS equations (see Appendix B of [8]), and attempts to solve them. The flavor D5- 
branes are taken to extend along the (x°, x 1 , x 2 , x 3 , if), r) directions and are smeared 
over the (9, <p, 9, (p) directions. These branes can be shown to preserve the same super- 
symmetry as the background for arbitrary values of the angles 9, (p, 9, (p [30]. Moreover, 
the smeared flavor branes will be sources for the RR 3-form, resulting in RR fluxes in 
the deformed background that can be observed as a "violation" of the original Bianchi 
identity. 
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The ansatz for the deformation of the MN background {H<$\ = 0, Ftg\ = 0) is, 

1 



ds 2 w = a'g s N c e^ 



a'g s N c 



dxl 3 + dr 2 + e 2h{r \dd 2 + sin 2 9dp 2 ) + 



e Mr) e 2fc(r) 

H — — ((wi + a(r)d9) 2 + (lo 2 - a(r) sm9d(p) 2 ) H ^— (^3 + cosOdipf 

The RR 3-form field strength reads 
a'g s N c 



(2.1) 



(3) 



4 



cui + b{r)d6) A (w 2 — 6(r) sin#<i(p) A (lu 3 + cos 9d(p) 



+b ; dr A {-d9 + sin 6dip A w 2 ) + (1 - K r ) 2 ) sin 9d9 A dp A co 3 ] 



(2.2) 



and automatically satisfies the Bianchi identity dF^ = 0. The left-invariant one-forms 
u n on S 3 are 



Ui = cos + sin ip sin l 
u 2 = — sin -^d^ + cos ip sin (9d<p, 
U3 = dip + cos ( 



We also introduce the new radial coordinate p, 

dp = e - k{r) dr, 

and the standard notation 

ei = d9, e 2 = sin 9dtp 
Ux — ojx + a(p)ei, u 2 = oj 2 - a(p)e 2 , 0J3 = UJ3 + cos 



(2.3) 



(2.4) 



(2.5) 



The metric then becomes 



ds w = e 



0(p)/2 



dxl 3 + e 2k ^dp 2 + e 2h ^( ei 2 + e 2 2 ) + 

e 2k(p) 



.-—(u 2 1 +u 2 2 ) + —(u 3 ) 



(2.6) 



where we have set a' = g s = 1 and N c has been absorbed into e 2h , e 29 , e 2k and dp 2 . 
The more familiar MN background is obtained from 2.1 and 2.2 with 

2r 



air 



sinh 2r 
e 2h(r) =rcoth2 r- 



b{r) = 0, k{r) = g{r) = 1 
r 2 1 



-2<f> 



sinh 2r 



2e~ 



2^0. 



Mr) 



sinh 2r 



(2.7) 
(2.8) 
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Note that even in the Nf = case, the system allows for other solutions -albeit with 
singular behavior in the IR- one of them is the solution with a(r) = often referred 
to as "abelian" 3 , others were studied in [28]. If we consider Nf ^ then there is a 
plethora of possible solutions. Among them there are solutions where the dilaton goes 
to a constant as r — > oo. Unlike MN that has a linear dilaton, the stabilized dilaton 
solutions do not correspond to a near horizon limit, gravity and field theory modes are 
coupled. However, as explained in [15] these are precisely the type of solutions we need 
to use as "seed" solutions for the rotation procedure to be used in Section 4, Appendix 
D. Note also that the ansatz for i*(3), (2.2), guarantees that Fcg\ is a closed form and 
thus it does not include any backreation yet. 

To include flavor branes and its backreaction, the action must be augmented by 
the DBI and Wess-Zumino actions for the flavor D5-branes. The complete action then 
reads 

S *Sgrav ~\~ ^sources; (2-9) 

where, in Einstein frame, we have 

Sgrav = / U ~ \{dM^) - 7^(3)) > ( 2 - 10 ) 

and 

N f 

Sources =T 5 V ( - / d & XC^ 2 ^9~6 + [ P[C (6) ] ) . (2.11) 
V J Me, JM 6 J 

where T 5 = and 2k\ = (27r) 7 . One of the effects of smearing the Nf — > 00 

flavor branes along the two transverse 2-spheres is that there will be no dependence 
on the angular coordinates (9,p>,9,p) of the functions f(p),h(p),g(p) and k(p) that 
determine our metric ansatz (2.6) — significantly simplifying the computations. After 
the smearing we can write 

S S ources = ^(-y d w x sin 9 sin ~9e* /2 + J Vol^(4) A C (6) J , (2.12) 

with the definition Vol(3^4) = sin 9 sin 9d9 A dip A d9 A dip. Once the smeared flavor 
D5-branes are incorporated into the background, the Bianchi identity for F(%) (which 
is identical to the EOM for C^\) gets modified to 

N f ~ Nf 

<iF(3) = —± sin 9 sin 9d9 A dp A d9 A dp) = A u 2 A e x A e 2 . (2.13) 

3 The abelian solution is characterized by a(r) = b(r) = 0, k{r) = g(r) = 1, e 2h ^ — r and 
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as a result of adding a Wess-Zumino term. This can be solved by adding the following 
term to the original F^) m (2-2), 



'(3) 



N f 

~—^ e i A e 2 A w 3 



(2.14) 



The full RR 3-form field strength now reads 



Nr. 



'(3) 



(tui + hex) A (u 2 - be 2 ) A tu 3 



+ b'dp A (-ei Awi + e 2 Au 2 )+(l-6 2 



AT, 



ei A e 2 A w 3 



(2.15) 



where the only modification is the appearance of the term involving Nf/N c in the 
second line. 

The first order BPS equations for this ansatz can be solved for b(p), h(p),g(p) leaving a 
system of three coupled differential equations for a(p), k(p) and <J)(p). The details can 
be found in [8, 28, 31]. 

In the present work we are interested in non-extremal generalizations of a particular 
type of flavored backgrounds: the ones with stabilized dilaton presented in [16]. The 
non-zero temperature breaks supersymmetry and thus we cannot resort to the BPS 
equations nor to the master equation formalism developed in [28]. We are bound to 
solve the second order equations of motion (EOMs) for b(p),h(p), g(p),a(p),k(p) and 
0(p) that will be obtained in the next section. 

2.2 Non-extremal flavored backgrounds 

One of our goals is to find non-extremal deformations of the flavored backgrounds 
presented in the previous section. To that aim we consider the following ansatz for the 
metric, 



ds 2 



IIB 

dsl 



- e * x dt + dx{ + dx 2 + dx 

J2k 



e 8x e 2k dp 2 

.29 



+ e^ 2 dsl 



(cu 3 + cos Odipf + e 2h {d9 2 + sin 



4 



+ — ((wi + a(p)d6) 2 + (oj 2 - a(p) sm6dip) 2 ) 



and the RR 3-form gauge field 



N r 



'(3) 



7(3) 



(2.16) 
(2.17) 
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where /( 3 ) is defined in (2.15). 

This ansatz is general enough to account for the effect of the backreaction of a large 
number Nf of flavor D5 branes, smeared along the 9, p, 8, (p directions, and spanning 
both the Minkowski and p directions. This gives a source contribution to the Fr$\ 
Bianchi identity, which we define as the smearing form 5(4): 

N f 

2(4)^^(3) = yWiAw 2 AeiAe 2 . (2.18) 

The factor e 8x ( p ' is a non-extremal deformation that accounts for the appearance of a 
horizon. The total Lagrangian for the bulk fields and smeared flavor branes is 

S 'S'grav ~\~ S sources , (2-19) 

where S'grav and Sources are given in (2.10) and (2.12) respectively. 

Before proceeding to present the EOMs note that the Wess-Zumino term in the 
flavor action does not involve the metric nor the dilaton so it will not enter Einstein's 
equations. Its effect is to change the equation for C( 6 ) which was d * = dF^ = 
and now in the presence of sources is dF( 3 \ = H( 4 ). The ansatz (2.16) assumes all 
the angular dependence is fixed by the symmetries of the background and the only 
non-trivial dependence is on the radial variable p . Following [25] we derive a one 
dimensional Lagrangian from which we can obtain the EOMs for all fields. We introduce 
the ansatze (2.16 - 2.17) in the action, integrate over the angular variables and drop 
the overall volume factor, 

Si = J{T - U)dr (2.20) 

where the explicit form of T and U and the details of the derivation are given in 
Appendix A. Note that this action should be supplemented with the constraint coming 
from reparametrization invariance, T + U — 0. 

Finally, defining s = -^f the EOMs read (seeAppendix A for details), 

8 4 v \jj 

+2 (g' + ti- 4x) 0' + 20' 2 + 0" = (2.21) 

2g'x' + 2h'x' - 8x' 2 + 2x'<p' + x" = (2.22) 
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± e -4( g+h -2 X ) e Ag + 16(} 4h _ e 4(g+k) _ ^(h+k) _ ^g g + g 2 

8 

+2e 29 (Ae 2h + Ae Ag+2h + e 2g+Ak - Ae 2h+Ak - e 2 °(-l + s)) a 2 - e As (-1 + e Ak ) a 4 

-Ae 2 %a (Ae 2h - e 29 (-l + s) + e 2 V) + b 2 (8e 2 ^ + Ae A9 a 2 ) 

+2k' (g' + ti - Ax' + 0') + k" = (2.23) 



^ e -2(g+2h) _ g e 2(g+h+k+Ax) + & 2g+Sx + e 2 9 +4fc+8rr + ^(h+k+ix) g _ 2 e 29+8:C S 

8 

+e 2 fl +8x s 2 + 2e 8x ^ e 2h + 2e 4 9+ 2ft _ ^fo+ft+k) _ e 2p+4fc + ^2^+4* _ + ^ Q 

+e W* (1 + e 4fc) fl 4 + 4e 8x 6 2 ^2h + e 2 3a 2^ _ ^Sx^ _ g 2 9( _ 1 + g) + ^2) 

+e 2^'2 + e 4 ff +2^/2 + 2/l , ^ _ 4x , + ^ + 2/l « + ^ = q (2>24) 



D -2(2g+h) 



32e 2{9+h+k+4x) + 16e 2h+8x + i6 e 2/l + 4fc + 8:r + 4 e 2 (9+ fc + 4a; ) s 



+4e 2^ & 2 _ ge 2 9+ 8x 6a _ 4e 2 9+ 8, ^ + ^ _ ^ fl 2 + ^,2 _ ^ 

+2</ (/>' - 4x' + 0') + 2</ 2 + g" = 



/2 



(2.25) 



-Ag-2h+Sx 



2e 2 °b 2 a + 6 (4e 2ft - e 2 %-l + s) + 3e 2 V) 



-a (4e 2/l (l + e 49 + e Ak ) - 8e 2 ^+ fe ) - e 2 ^ +4fc - e 29 (-l + s) + e 2g (l + e Ak ) a 2 ) 
+2a' (2g' - Ax' + 0') + a" = (2.26) 



-Ae 8x b + Ae 8x a + e 29 ~ 2h+8x a (l - s - 2ba + a 2 ) + b' (-8a;' + 20') + 6" = 0. 



(2.27) 



These equations will be solved together with the constraint coming from reparametriza- 
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tion mvariance, 



]_ e -4(g+h) 



a 



256 

_ e <ig+8x _ lQ e <±h+8x _ g e 2(g+h+k+4x) s + 2 e 45l+8:E S - e 49+8a: S 2 

_ 2e 2g+8 : c ^2/x _j_ 4e 4 9 +2/i _ ^(g+h+k) _ e 2g+Ak + ^h+Ak _ g 2g^_ 1 + g ^ „2 
_ e 4 9 +8x ^ + e 4^ fl 4 + 4e 2 9 +8x a ^2*. _ & 2 g ^ + ^ + ft 

-4e 29+8:E (2e 2/l + e 2 V) 6 2 + 2e 69+2 V 2 + 2e 2(s+/l) 6 /2 

-i (g' 2 + h' 2 + 2/i / (fc / - Ax') + 2#'(2/i' + Jfe' - 4x') - 4JfcV) 
+l(-2g' - 2ti — k' + Ax')(j)' - \<j>' 2 = 0. 



(2.28) 

3 New non-extremal flavored solutions with stabilized dilaton 

In this section, we numerically solve the Einstein equations of motion (2.21)-(2.28). 
Our method combines the virtues of the approaches developed in [23] and [24, 32]. 
Following [23] we first find a set of parameters that completely determines the UV 
behavior of the solutions to any order. In our case there turn out to be eleven free 
parameters. The IR behavior is determined by regularity at the horizon and involves 
seven free parameters. We then solve numerically the equations of motion and the 
constraint using these expansions as boundary conditions. To derive the EOMs we 
follow the framework in [25] (see Appendix A for details). 
Two points are worth mentioning 

• The solutions we are looking for are non-extremal deformations of the supersym- 
metric flavored solutions studied in [16]. Therefore, we fix some UV parameters 
so that the UV asymptotics reduce to the supersymmetric solutions when the 
non-extremal parmeter goes to zero. 

• An important feature of our solutions is a dilaton whose value is "stabilized" - 
i.e. it asymptotes to a constant value in the UV. This property is required if the 
rotation procedure is to produce new solutions where the gravity modes decouple. 
Recall that, as discussed in [15], [16], in the Nf — 0, T = case we have a two- 
parameter family of solutions. MN corresponds to the particular case when both 
parameters are taken to zero -or equivalently when the decoupling limit is taken- 
while we are interested in solutions where one of the parameters is infinity and 
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the other one is finite. Thus in the Nf — y 0, T — > limit our solutions are not 
expected to reduce to MN. 4 



3.1 UV expansions 

To determine the UV behavior of the metric and gauge functions for a finite-temperature 
deformation of the flavored solutions of [16], we begin by taking a UV series expansion 
of the form 

oo i oo i 

e 2h =J2Yl h ^ f 9 eH =HH fid f 9 e 4(1 ~ iW3 

i=0 j=Q i=l j=Q 

oo i oo i 

oo i oo i 

i=0 j=0 i=l j=0 

oo i 

i=l j=0 

and requiring that it satisfies the Einstein equations given in Appendix A. For an 
expansion of this form, we find that all of the coefficients — up to arbitrary order — are 
determined in terms of a set of eleven free parameters. Details of the relations between 
the expansion coefficients are given in Appendix B.l. We next set a number of the free 
parameters to agree with the supersymmetric UV asymptotics of [16] (summarized in 
Appendix B.2): 



/i,o = 1, / 3 ,o = 3/4^ - 3s/Ac 2 + + 51s 2 /84, hi, = -1/4 + 13s/32, 

, 0/0/ -512 C 3+8c_+3s(32-32s+17s 2 ) , , , . 

fc ,o = 2c+/3, fc 3i0 = ± fti,i = 1/2 - s/4, 

a 2 , = 0, a 4 ,o = 2, 64,0 = (3.2) 

After doing this, the asymptotics reduce to those of [16] in the limit where the 
parameter a^o goes to zero. In order to agree with the convention used in [21], we will 
refer to the parameter £5 as C2 in what follows. 5 With the choice of parameters in 

4 Using the notation of [16] let (3 and c be the parameters that characterize the family of solutions. 
We have the following cases: l)when (3 — c — we get MN 2) when /3 — » 00, c->oowe obtain KS, 
3)when j3 — > 00 and c is a free but finite parameter we get the KS baryonic branch 

5 Note however, that the sign of our C 2 is flipped with respect to theirs: we expand e 8x , whereas 
they expand e~ 8x . 



-12- 



equation (3.2), the UV asymptotics of our metric and gauge functions are given by 

e 8 * = 1 + C 2 e- 8 ^ - + (e- 8 ^ 3 ) 

2c_|_ 

p 2k _ 2 p4p/3 ,s_ e~^ 3 (400 - 400s + 91s 2 - 160(-2 + sfp + 64(-2 + sfp 2 ) 8p/3 
~3 + + 2 96c + + 1 J 

e 2h = \c + e^ 3 + l(-8 + 13, - 8(-2 + s)p) 

+ -^e~ 4 " /3 (208 - 208s + 43s 2 - 64(-2 + sfp + 64(-2 + s) 2 p 2 ) + £>( e - 8p/3 ) 
256c+ v ' ' 

e 2 ^ = c+e 4 ^ 3 + (l + ^ + (-2 + s)p 

+ e~ 4p/3 (208 - 208s + 43s 2 - 64(-2 + s) 2 p + 64(-2 + s) 2 p 2 ) + C(e" 8p/3 ) 
64c + v ' 7 

e 4 = l _ 3e~ 4 ^ 3 3 + e -8p/3 ( 6 + s 2 (51 - Up) - 48p + 6s(-l + 8p)) + G(e _ 4p) 

c_|_ 8c^_ 

a = 2e- 2 " + O(e- 10p/3 ) b = 2e~ 2f> (2 - s)p + 0( e - up/3 ). (3.3) 

All of these expressions in fact depend on C 2 , but only at higher order than we have 
shown here. The parameter c_ also enters at higher order. In our numerical solutions, 
we will set c_ equal to for simplicity. This choice of UV behavior leaves us with three 
free parameters: c+,C2,and s. 

3.2 IR asymptotics 

Near the horizon, p = p h , we expand the equations of motion in a power series up to 
fifth order: 6 

e -8x( P ) = Xl ( p _ Ph ) + X2 ( p _ Ph f + X3 ( p _ Ph f + X4 ( p _ Ph y + X5 (p _ Ph f + . . . 

e 2h ^ = h + h 1 (p- p h ) + h 2 (p - p h ) 2 + h 3 (p - Phf + h 4 (p - p h ) 4 + h 5 (p - p h ) 5 + ■■■ 

e Am = f + A( P - Ph ) + f 2 {p - phf + Hp - Ph f + Hp - Ph f + Hp - Ph f + ■■■ 

a(p) = a + ai(p - p h ) + a 2 (p - p h f + a 3 (p - p h f + a A (p - p h f + a 5 (p - p h f H 

b(p) = b + hip - p h ) + b 2 (p - p h f + b 3 (p - Phf + h(p - phf + b 5 (p - phf + ■■■ 

(3.4) 



6 Thc horizon expansion could be taken to higher order to improve accuracy. However, due to 
the presence of s the expressions grow unwieldy in complexity. We find that a fifth order expansion 
provides enough accuracy to obtain well behaved numerical solutions. 
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Demanding that these expansions satisfy the equations of motion, we derive expressions 
for x 2 ...,hi...,gi...,ki...,fi...,ai...,b 1 ... in terms of x 1 ,h ,go,k ,ao,b Q . 7 Along with / , 
this gives seven independent parameters coming from the horizon expansion. The 
expressions for the dependent coefficients are quite cumbersome and we will not present 
their explicit form here; a Mathematica file is available from the authors upon request. 

3.3 Numerics 

In this section we outline the numerical procedure used to find solutions to the equa- 
tions of motion (2.21)-(2.27) that have a horizon at finite p = ph and obey boundary 
conditions (3.3). Due to the large dimensional parameter space (3 free parameters in 
the UV and 7 in the IR) this is a difficult numerical problem. The numerical method 
we use does not pretend to solve the problem in all generality and is not a tool to 
thoroughly study the full parameter space. However, within certain limitations to be 
elucidated below, it allows us to construct numerical solutions with enough accuracy 
to determine the temperature, energy and specific heat of these backgrounds. More 
details of the numerics are given in Appendix C. 

Our strategy is to start in the UV where there are fewer free parameters. We pick 
a set of values for c + and C 2 , 8 which determine the UV behavior of the background. 

With these UV boundary conditions, we integrate back from p^ and look for values 
of c + and C 2 which produce solutions with horizon behavior. We will refer to these 
solutions as the UV-shot solutions. The value of p^ is chosen such that the dilaton 
reaches its asymptotic value, e A ^\ Pac ~ 1. We next require that this numerical solution 
matches the horizon expansion (3.4) and its derivatives up to fifth order. Doing this 
for a given choice of (c+, C 2 ) determines the free parameters (xi, h , g , k , f , a , bo). A 
black hole solution will not exist for every choice of (c+, C 2 ); for some values of (c, C 2 ) 
there will not be a horizon and for others there may be naked singularities outside the 
horizon. 

In order to match the solution to the horizon expansion, we define a "mismatch" 
function evaluated at a point po close to the horizon, p = Ph + e. 

m(po) =[x sh (p ) - X num (p )] 2 + [x' sh (p ) - x'numiPo)} 2 + [fi^(Po) ~ SW(po)] 2 + 
Wsh(Po) - g'num(Po)} 2 + [MPo) ~ K urn (p )} 2 + [h' sh (p Q ) ~ U num (p )] 2 + 

[ksh(Po) - k num (p )} 2 + [k' sh (p ) - k' num (p )} 2 (3.5) 

The subscript num denotes the UV-shot numerical solution and sh denotes the series 
expansion (3.4). In [21] the authors used a similar technique to find abelian, a = b = 

7 The higher order coefficients do not depend on /o. This corresponds to the invariance of the EOMs 
under constant shifts in the dilaton. 

8 For simplicity, we set c_ equal to zero. 
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0, solutions. Unlike [21], here we are dealing with non-abelian solutions and due to 
numerical inaccuracies near the horizon, including a, b, and in the mismatch function 
does not allow for an accurate matching. Our strategy to find the horizon parameters 
will be a three step procedure. First we will minimize a mismatch function m (that 
does not contain a, b, and <p) and determine (xi, h , g , k ). Then we will tune a$ 
and bo to eliminate unwanted behavior in b. Finally, we will use the invariance of the 
EOMs under constant shifts of the dilaton to set fo so that the horizon-shot and UV- 
shot dilaton agree at the UV boundary. Keep in mind that we want to perform the 
matching as close to the horizon as possible and obtain stable solutions. We use the 
word "stable" loosely to mean horizon shot solutions that agree with the UV shot ones. 
This is not an easy task because close to the horizon the UV-shot solutions always 
display numerical inaccuracies. To overcome this problem we do the following: we first 
match at a larger value of e = 0.7, then we use the resulting values of the matched 
horizon parameters as seeds for a new match performed at a lower e = .15. In the 
last step we also constrain the allowed variance of (xi, g , h , k ) around the seed values 
until the resulting e = .15 match gives horizon shot solutions which agree with the UV 
shot solutions. The end result is m(p/ l + .15) < 1CT 7 and a set of parameters highly 
tuned to give stable solutions. 

With these parameters we shoot from the horizon and show that our solutions also 
satisfy the constraint coming from reparametrization invariance, T + U = 0. We use 
WorkingPrecision = 40 and obtain solutions that satisfy T + U < ~ 10~ 6 throughout 
the interval. Figs.(l - 4) show two solutions for different (c+,Cz) values. 

Let us comment on the relation between the numerical solutions presented here 
and some known cases. The non-extremal non-flavored with stabilized dilaton exists 
in the literature only for a(p) = [21]. We can reproduce the results of [21] with 
improved numerics. The constraint is much better satisfied now, T + U < 10~ 6 here as 
opposed to T + U < 10~ 3 in [21]. In addition, the mismatch function is of order 10~ 7 
here and was 10~ 4 in [21]. The extremal flavored solutions with stabilized dilaton [16] 
used a slightly different framework, they reduce the problem to an unknown master 
function 9 for which they solve numerically. This makes difficult a direct comparison 
of the numerics. However, we can compare with their UV and IR expansions for 
e 29 , e 2h , e 2h , etc. Since we use their UV deformed by a non-extremality parameter C*2, 
our UVs automatically reduce to the ones in [16] when = 0. The IR is a bit more 
delicate. Their solutions are singular in the IR, e 29 , e 2h , e 2k , go like logr as r — > 0. 
If we set C 2 = in our code we run into a singularity at p — > and reproduce the 

9 The master function is obtained form the BPS equations and thus, is not a formalism we can use 
in the non-extremal case 
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solutions of [16]. Also note that we find regular horizons only for very large values of 
the non-extremality parameter; for small values of C2 we run into singularities. 




lih " G 7 8 " 9 rh 5.2 5.4 5.6 



Figure 1. Horizon-shot and UV-shot solu- Figure 2. Horizon-shot and UV-shot solu- 
tions for metric functions at s = 1, c+ = 3, tions for metric functions at s = 1, c+ = 3, 
C*2 = 800000. e k - orange, e g - blue, e h - pur- C2 = 800000. Zoomed into region near hori- 
ple, e Sx - pink zon at ph ~ 5.097. 




Figure 3. Horizon-shot and UV-shot solu- 



tions for metric functions at s 
C 2 = 5000. 



1, c+ 



Figure 4. Horizon-shot and UV-shot solu- 



50, 



tions for metric functions at s = 1, Cj 



50, 



Ci = 5000. Zoomed into region near horizon 
at ph ~ 3.19 



3.4 Temperature of the solutions 

The expression for the horizon temperature of our solutions follows from the standard 
prescription of analytically continuing the solution to Euclidean time and imposing the 
periodicity of this coordinate. The temperature is then given by T = 1/(3, where the 
period (3 is determined by requiring the regularity of the Euclideanized metric at the 
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horizon. We use the numerically determined parameters described in section 3.2 to 
describe the metric at the horizon; of these parameters, only x\ and ko enter into the 
expression for the temperature, which is found to be 



In [21] it was shown that the horizon temperature is the same for pre- and post-rotated 
solutions. It can be checked that this is still the case here. 

In Figs. (5) and (6) we see that, for fixed values of c + and C 2 , the horizon tem- 
perature has negligible dependence on s e Nf/N c . As we vary < s <~ 10 the 
temperature remains constant; we conclude that the temperature is not affected by the 
backreaction of the flavor branes. 

In Figs. (8), (7) and (9) we show how the horizon temperature varies with C 2 for 
fixed values of c + and s. In Fig. (10) we present the temperature as a function of the 
horizon area. Note that as G% decreases the temperature increases indicating that it 
diverges in the extremal limit. Also, for C 2 — » oo figure (9) shows that the temperature 
goes to a very small but non-zero value, T c . This behavior is reminiscent of the wrapped 
five branes black holes of [25]. In that framework T c is the Hagedorn temperature of 
the little string theory and is also a critical temperature at which a first order phase 
transition occurs (chiral symmetry restoration and deconfinement) [25]. A similar phase 
transition could exist here. If that is the case, the U-dualities described in the following 
section will probably map it to a phase transition in Klebanov-Strassler. This is an 
issue that deserves further study. 
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Figure 5. Temperature at the horizon versus Figure 6. Temperature at the horizon versus 
N f /N c . c+ = 50, C 2 = 5000. N f /N c . c+ = 3, C 2 = 800000. 




Figure 7. Temperature at the horizon versus Figure 8. Temperature at the horizon versus 
C2. c + = 50, s = 1. C2. CLj_ = 3, s = 1. 
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Figure 9. Temperature at the horizon versus Figure 10. Temperature at the horizon ver- 
for very large values of C2. c+ = 3, s = 1. sus -A/V^. c+ = 3, s = 1. 



4 Non-extremal backgrounds with flavored resolved deformed 
conifold asymptotics 

In [15] the authors presented a solution generating technique that takes a IIB back- 
ground of .D5 branes wrapped on the S 2 of the resolved conifold (non-trivial and 
0) and produces a background wich also has D3 charge (non trivial Ftg\, F/g\, and 
0). After a scaling limit the resulting geometry represents the baryonic branch of the 
Klebanov-Strassler solution. The initial space has topology R 1 ' 3 x A4 e and preserves 
M = 1 supersymmetry . The algorithm consists of performing three T-dualities in the 
R 3 directions, lifting to M-theory, boosting with rapidity (5 in the eleventh direction, 
reducing back to ten dimensions and finally T-dualizing back in R? . This procedure is 
equivalent to a rotation in the space of Killing spinors. We will therefore refer to this 



algorithm indistinctively as a chain of dualities or as a rotation. 

In the previous section we found non-extremal solutions representing a background 
of N c wrapped color D5 branes and Nj smeared D5 flavor branes (Nf/N c ~ 1). In this 
section we want to apply the procedure developed in [15] to these backgrounds. 

After applying the rotation procedure, taking /3 — > oo to decouple the gravitational 
modes, and performing the rescalings as described in [21], we obtain the transformed 
solutions in Einstein frame (see Appendix D for details): 



dsj IB = N c 



e -0/2 H -l/ 2( _ e -8x rft 2 + dx , dx i^ + e 30/2^1/2 rfs 2 
4 'V J )' ^(3)-- e -4 



*<3)=~/(s), ^(3) = -e- 4 ^e^* 6/(3) 



F (fi) = -Ni{l + * 10 ) 



Vol(4) A d 



e -Ax 



(4.1) 



where 



H 1 ' 2 = Ve~ 2 * - e~ to (4.2) 

and the explicit forms of /(3) and dsf are given in equations (2.15) and (2.16) re- 
spectively. Note that, unlike more standard finite temperature solutions where the 
non-extremality factor e~ 8x enters only in the g rr and gu elements of the metric, here 
the warp factor H and the fluxes and depend on e~ 8x . 10 

The use of U-dualities is a well stablished solution generating technique in super- 
gravity. In the presence of smeared flavored branes one could question the validity and 
meaning of the uplift to eleven dimensions. Two pieces of knowledge provide useful 
insight on this issue: 

1. In [33] the authors studied the uplift of smeared D6 flavor branes to M-theory. 
They argued that the violation of the Bianchi identity in ten dimensions implies 
that the eleven dimensional-geometry will no longer be Ricci flat; it is no longer 
of G2 holonomy but it does carry G 2 structure. Thus, the flavors appear in 
eleven dimensions as intrinsic torsion. This eleven dimensional theory is no longer 
maximally supersymmetric (it is 1/8 BPS) and therefore it is no longer unique. 
This is an interesting possibility for the interpretation of the uplift of smeared 
flavor branes. For us, however, the important point is that once the reduction to 
ten dimensions is carried out the correct smearing form is recovered. 



10 Note also that in [21] the factor of e~ 8x in H/$\ trivially cancelled with a similar one coming from 
the six dimensional Hodge dual. This is not the case here. The reason is that in [21] the authors 
considered a particular point in the baryonic branch (a = 0) while in the present work we consider 
the general case, a(p) (see equation (D.4)) 
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2. In [16] the authors studied a transformation to generate SU (3) structure solutions 
of IIB supergravity starting from non-Kahler backgrounds describing wrapped 
D5 branes with additional flavor (smeared) D5 branes. Working with the BPS 
equations and the ST/ (3) structure of the backgrounds they presented a solution 
generating technique that amounts to a rotation in the space of Killing spinors. 
This rotation procedure is well defined (even in the presence of smeared flavor 
branes) and is equivalent to the chain of U-dualities alluded to above. 

In the present work we are dealing with a non-extremal deformation of [16]. How- 
ever, since the backgrounds we are working with are not supersymmetric, we cannot 
use the BPS formalism of [16]. On the other hand, we are not interested in the inter- 
pretation of the background in eleven dimensions, so we will not attempt to address 
the questions presented in [33]. Instead, guided by the results of [16], we will forge 
ahead and apply the chain of dualities to the flavored, non extremal backgrounds with 
stabilized dilaton found in section 3. However in order to claim that the outcome of 
these dualities is still a solution of IIB plus sources, we have to show that the rotated 
backgrounds are indeed a solution of the EOMs. 

4.1 Equations of motion with D3 and D5 charges and smeared flavor branes 

In this section we will verify that given the EOMs before the rotation procedure is 
applied, the EOMs after the rotation are also satisfied. The EOMs before rotation 
were presented in section 2.2 and derived in Appendix A. After the rotation procedure 
is carried out, the supergravity background contains if( 3 ) and F( 5 ) fluxes in addition to 
the _F(3) that was present before rotation. The full type IIB action, supplemented by 
the DBI action for the smeared flavor branes, is then given by 



S — Sjjb + Screes (4.3) 



where 



Sub = ttV / ^9 R - / ( dtp A *d<p + e-^3) A *H {3) 

ZK 1Q J 4K 10 



and 



+ e<%) A *F (3) + ^F (5) A *F (5) - C (4) A F {3) A H {3) j 
Sources = J ( Vol (4) A dp A dip e^ /2 ^\g ab + e-^ 2 B ab \ 



(4.4) 



- C(6) + C(4) A B(2) ) A =(4) 



(4.5) 
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Note that the metric is in Einstein frame and that g a b and B a ^ respectively represent 
the pullbacks of the metric and B^) to the flavor brane world volume. Recall that 
before rotation = and thus, the DBI action was only the pullback of the metric. 
After rotation we need to specify the form of the NS potential, -8(2)- in the extremal 
case (T = 0), this potential should reduce to the one in [16]. We propose the following 
ansatz: 

B{2) = bx(p)<jJ 3 A dp + b 2 (p)e 1 A e 2 + b z e x Aw 2 + b A (p)e 2 A u)i + b 5 (p)u l A Q 2 . (4.6) 

To make the gauge degrees of freedom more apparent it is convenient to parameterize 
5(2) in a slightly different way. Let us introduce b 2 (p) = b 2 (p) + (1 — a(p) 2 )b 5 (p) . It is 
then easy to verify that 



and that b 5 is just a gauge degree of freedom. Demanding that H^) — d[B^}, we get 
5 equations. This system of equations can be reduced to just one differential equation. 
We choose a gauge such that &i(p) coincides with the one in [16] (see Appendix E for 
details). 

The flavor brane world volume is dt A dx\ A dx 2 A dx$ A dp A dip = Vol(4) A dp A dip, 
and the smearing form is given by 5(4) = | sin 9 sin 9 d6 A dip A d6 A dip, where we have 
set a' = g s = 1. As was the case before rotation, the C(e) A 5(4) term in the source 
action vanishes. Note that in the EOMs that follow from this action, the contribution 
from the Wess-Zumino term in the source action is exactly cancelled by that from the 
Chern-Simons term in the bulk. 

The Bianchi identities for the RR and NS-NS forms after rotation are 



The first is unchanged from before the rotation, and the latter two can be checked to 
hold given the EOMs before rotation. Note that this does not give us the full expression 
for -6(2), but only the components which are orthogonal to 5a), namely B p ^. However 
-8(2) only appears in the source action, and there it is either wedged with 5(4), or pulled 
back to the flavorbrane world volume (which is orthogonal to S( 4 )). Therefore only this 
orthogonal component will be relevant for our equations of motion. A more detailed 
derivation of the general form of B^ is given in Appendix E. 
The EOMs for the fluxes after rotation are 




dFQ) = 5(4), 



dH m = 0, 



dF^ + F( 3 ) A if( 3 ) = -B( 2 ) A 5(4) 



(4.8) 



d (<T* * H {3) ) = F (3 ) A F (5 ) + 



H 



Vol( 4 ) A 5(4) 



(4.9) 
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d(e**F (3 )) =-if( 3 )AF (5) (4.10) 

where the term with T-L in the first equation comes from varying the source term. These 
can be checked to hold given the EOMs before rotation. 
The dilaton EOM is 



d (*#) = y F (3) A *F (3) - ^"#(3) A *if (3) + (2k? ) . 



where 



,0/2 



50 



\9ab + e~^ 2 B ab \ 



1 5(4) | e 



V\9ab + e-WBab] 



Vol (4) A dp A dip A S(4) 
(4.11) 



and can be checked to hold given the EOMs before rotation. 
Lastly, we have the Einstein equations 
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(4.12) 



where the variation of the Chern-Simons term in the bulk action is cancelled by the 
variation of the Wess-Zumino term in the source action. Here, T*° UTCes is obtained from 
(after rotation) 
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(4.13) 
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(4.14) 



Again, the Einstein equations can be checked to hold given the EOMs before rotation. 
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4.2 Rotating the solutions 

It is now a simple procedure to use equations (4.1) along with our numerical solutions 
to produce numerical solutions for the background after rotation. Before presenting 
and commenting on the numerical results let look at the large p behavior. 

4.2.1 Asymptotics after the rotation 

In order to write down the form of the metric the UV after the rotation, we identify 
from (4.1), 



Qxx = e " H 1 g t t = e ft ' g pp 



e 29 „ / ,,, / ,e a 

T 



(4.15) 



Using (3.3) and (4.15) we obtain, 

^16c 2 C 2 +3 (-2 + s(2 - 16^) + 16p + s 2 (l + 4p)) J 



^e 2 "/ 3 e -2 P /3 
7s 2VScTs" 3/2 

^e 2 ^ 3 p -2p/3 



\/— L e l/ — 2p/3 / \ 

-S« - 24^.3/2 (l6c 2 + C 2 +3 (-2 + s(2 - 16p) + 16p + s 2 (l + 4p)) J +■ ■ 

9 PP ~ y^e 2 ^ 3 v ^+ 24 ^_^ / - ^-6 + 164C 2 + g (6-48p)+48p+3 g 2 (-ll+4p)j +■ • • 

e -2p/3 / \ 

i v ^+— =^ -6 + 164C 2 + s 2 (9-36p)+48p + 6s(-7+8p) +••• 
64 A /6c+ v / s V / 



In p-2p/3 
^+ e 2p/ 3i - ' 

2 

Q c g-2p/3 / \ 

— ± e 2p/3 y /s+ rA = ■= -6+16c 2 C 2 +48p-18s(-3+8p)+s 2 (-39+60p) +• • • 
2 64y/6c + y/s \ J 

~ V ^ + 96 ^-^- ^-6 + 164C 2 + g (6-48p)+48p + 3 g 2 (-ll + 4p)j +■ ■ ■ 

(4.16) 

It is clear that this expansion is not valid as s — > 0. However, this does not mean 
that we cannot obtain the expected KS asymptotics in that limit. It just indicates 
that we should take the s — >■ limit before performing the rotation. Indeed, setting 
s = before the rotation in (3.3) and defining c+y/8A(p) = a/24 P + 8C2C+ — 3 to avoid 
clutter we obtain the following metric after rotation 11 , 

11 The effect of flavor in the UV can be heuristically understood by looking at the warp factor 
H 2 = e- 2 ^ - e~ 8x and at the expansions (3.3). In the flavored case H 2 ~ e" 4 ^ 3 ^- + e>(e~ 8p / 3 ) while 
in the unflavored case V 2 ~ p(C2, c+)e~ 8p / 3 + 0(e~ 4p ), where p(p; C2, c+) denotes a polynomial in p. 
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4p/3 e -4p/3 / 

+ 3 (-847 + 3504p - 3456p 2 + 2048p 3 ) ^ + C(e" 8p/3 ) 



4p/3 e -4p/3 / 

+ 3 (-847 + 3504p - 3456p 2 + 2048p 3 ) J + C(e" 8p/3 ) 
9 PP ~ -fA(p) + 0(e- 8 ^) g ee ~ -±A(p) + —^A(p)(2p - 1) + O^) 

1 r p -4p/3 

^ ~ -c + A(p) + 0(e- 8 ^ 3 ) ^- ~ -±A(p) - _,l(p)(2p - 1) + 0(e- 8 ^ 3 ). 

(4.17) 

Using a variable u = e 2p//3 , the metric reads at leading order 

ds 2 = -^-(d Xl dx l ) + + ds 2 T11 ) + ■■■ 0(u- 2 ), (4.18) 



where A(u) = y log it + fC^c^ — ^ contains the typical log behavior of KS asymp- 
totics. . Note that, since numerically the work of finding non-extremal solutions is done 
before the rotation it is trivial for us to set s to any value, however small, and then 
rotate. 

4.2.2 Rotated flavored solutions 

After rotating the numerical solutions obtained in Section 3 we obtain Figs. (11 - 18) 12 . 
Note that, 

• The addition of flavor alters the behavior of the dilaton and the g XX i9tti 9pp parts 
of the metric in the UV, but not in the IR. On the other hand, flavor affects the 
metric of the compact manifold at all energy scales, as shown in Figs. (17, 18). 

12 The UV behavior of the rotated solutions is extremely sensitive to the precise UV behavior of the 
seed solutions, due to the H factor. For this reason, we have generated these plots from the UV-shot 
seed solutions. 
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Figure 11. g pp metric element for solutions Figure 12. g pp metric element for solutions 
after rotation, c+ = 50, C2 = 5000. From after rotation, c_|_ = 3, C2 = 800000. From 
bottom to top: s = 0, 1, 19/10, 51/10, 91/10. bottom to top: s = 0, 1, 19/10, 51/10, 91/10. 




Figure 13. g xx (red) and gu (blue) metric Figure 14. g xx (red) and gu (blue) metric 
elements for the unflavored solution after ro- elements for the unflavored solution after ro- 
tation: c+ = 50, C 2 = 5000, s = 0. tation: c+ = 3, C 2 = 800000, s = 0. 



• Comparing g 9e and g§g for s = 61/10 and s = 2/5 in Figs. (17, 18) we can see 
that the behavior of the curves has been swaped. For s = 61/10, grows as we 
approach the horizon while gee decreases. For s = 2/5 it is the other way around. 
This swapping occurs as we cross s = 2. 

• From (4.1) we can see that after the rotation the compact part of the metric is 
still an expanding T 1,1 . This is also shown in(17, 18) where we can see that for 
large p we have gee ~ 9ee ~ fflW- 
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Figure 15. g xx (red) and gu (blue) metric Figure 16. g xx (red) and g tt (blue) metric 
elements for flavored solutions after rotation: elements for flavored solutions after rotation: 
c + = 50, C 2 = 5000. From top to bottom: c+ = 3, C 2 = 800000. From top to bottom: 
s = 0, 1, 19/10, 51/10, 91/10. s = 0, 1, 19/10, 51/10, 91/10. 




SSR / 
SSIV-j 



Figure 17. \g^ (blue), gee (red), and g§g (yellow) metric elements for 
flavored solution after rotation. The different groups correspond toi s = 
(solid), s = 2/5 (dotted) and s = 61/10 (dashed). c+ = 3, C 2 = 800000. 

5 Energy and specific heat 

In order to assess the thermodynamic stability of our solutions, we need to determine 
the specific heat, C v . This is obtained from the expression for the ADM energy of our 
solutions and their temperature via the standard thermodynamic relation C v = dE/dT. 
T is the numerically determined temperature given in equation (3.6), and E is given 



- 26 - 



UIR 



Figure 18. (blue), gee (red), and gss (yellow) metric elements for 

flavored solution after rotation. The different groups correspond to s = 
(solid), s = 2/5 (dotted) and s = 61/10 (dashed). c+ = 50, C 2 = 50000. 

by the expression for the conserved ADM energy [34] (where 167rGio = 2/^ = (27r) 7 ) 

Eadm = -^4r- / V\9^\( S K - 8 K Q )dS?. (5.1) 

Here, 8 if and 8 i^o are the extrinsic curvatures of two 8-dimensional submanifolds St of 
St at constant p. £f is a constant-time 9d spatial slice of the entire geometry, and we 
take for S t the boundary manifold as p — > oo. S K corresponds to our finite temper- 
ature backgrounds, whereas 8 K is the extrinsic curvature of a reference background 
[34, 35]. For the reference background, we take one of the family of zero-temperature 
BPS solutions discussed in [16, 28]. The UV asymptotics of these backgrounds are 
summarized in Appendix B.2. 

Our Euclideanized lOd metric reads 

ds 2 = goodr 2 + g xx dxidx l + g pp dp 2 + gee{e\ + e\) + g§§(ul + cl>|) + 5y>V> (5.2) 
The metric on a constant time slice S t , and on its 8-dimensional boundary S^° are 

ds| t = g xx dxidx l + g pp dp 2 + g e e{e\ + el) + g§§(uj 2 + toj) + g^ (5.3) 
ds 2 sr = g xx dxidx l + g w {e\ + ej) + g§g(uif + G%) + g^ u\, (5.4) 



where the metric functions in Einstein frame are 
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before rotation, and 
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_ 1 p 30 s /2o/l/2 2fc fl 



(5.6) 

after rotation; l B' subscripts indicate the BPS backgrounds. Here H = (e~ 2( ^—e~ 8x ) and 
U B = (e~ 2 * fl -Ki{), where k>2 is a free parameter with < k 2 < e that parametrizes 
the rotation, see [16]. In terms of the lid boost parameter (3, we have k 2 = e~*°° tanh 0. 
Note that for the BPS solutions e~ 8x = 1, and we have also allowed for an arbitrary 
rescaling of the r coordinate via a constant parameter ub, reflecting the freedom in 
choosing the period /3 of the Euclidean time coordinate for these backgrounds. 

It is now straightforward to calculate the extrinsic curvature of an 8-manifold St 
at constant p. We obtain 

8 K bef = e -^-KP)-MP) (2 0'(p) + 2 g\p) + 2h'(p) + k\p)) (5.7) 
for the background before rotation, and 

* K ft - 1 VM> r -'^-k(p)-4x(p) 

aft 2H(p) 5 / 4 

+ n^ e ~^~ Hp) ~ Mp) (w'(p) + zg'(p) + 2h'(p) + fc'(p)) (5.8) 
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for the rotated background. 13 
The expressions for dSt are 
1 



dS-, 



tbef 



^ e 2 9+ 2h+k dv A ds 



Me 



dS t 



of I 



_ n l/2 e 3<f> e 2 a +2h+k dV3 A (5>g) 



where dV% = dx\ A d%2 A dx^ is the volume form on the Minkowski spatial directions, 
and dSMe — (sm9sin9)d9AdpAd9AdipAdip is the volume form on the compact cycles. 
Putting all of these together and integrating over the compact 9, (p, 9, <p, ip directions 
we obtain expressions for the ADM energy before and after the rotation: 
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(5.11) 



p— >oo 



Here V3 = Vol(-R 3 ), and we have expressed the results as energy densities per unit flat 
3- volume. 

We evaluate these expressions at a large but finite value of p = p* by inserting the 
UV asymptotics of the BPS and finite temperature solutions. The BPS asymptotics are 
given in terms of the free variables Q Q , c+, c_, po and /^o, and are listed in Appendix B.2. 
The finite temperature asymptotics were given in section 3.1, and are a deformation of 
a BPS background with Q Q = —N c + Nf/2 and po = 0. 

In order to calculate the ADM energy, we first match the geometries and matter 
fields of the BPS and finite temperature backgrounds at p*, this amounts to match- 
ing the metric functions goo, 9xx, gee, g§§, g^, as well as the dilaton. Only after this 
matching is performed we let p* — > 00. 

We begin with the case before rotation. To avoid confusion with the parameters c + 
and c_ in the finite temperature solutions, we will write the BPS c + and c_ parameters 
in bold as c + and c . First we pick some fixed values of c + ,C 2 and s to specify the 
finite temperature background. We begin the matching by using the freedom in v B to 



set v B = e 



Next we match e^ 2 = e* s/2 to 0(e~ 12p * /3 ) by setting 



-12p»/3 



C2S 



(5.12) 



13 



The BPS backgrounds will have e Ax = 1. 
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As a result, goo an d g X x are now matched to 0(e 12 ^*/ 3 ). Now, we match e 2k = 
e 2kB ,e 2h = e 2hB , and e 2g = e 2gB , respectively, to 0(e~ 8p * /3 ), by setting 

c + = c + + e- 12 ^ 1 -^ 
+ + 48 

Qo = -1 + J - 4 e ~ 8p * /3c, 2(-5 + 4p*)(-2 + s) 
2 lb 

c_ = ^4(-llC 2 s + 4C 2 sp*), (5.13) 

which results in the matching of g$g, g^, g^ to C(e- 8 ^/ 3 ). 

With these definitions, the dilaton and all of the metric functions goo, 9d, gee,ggg, g^ip 
are matched up to (at least) (9(e~ 8p */ 3 ). The freedom in the po parameter in the BPS 
solutions is not needed, and we set it equal to zero as was done for the finite temperature 
solutions. Taking p* — > oo, the resulting ADM energy density is finite, independent of 
s, and equal to 

e bef = 5 -^ (5.14) 
; 96vr 4 v ' 

After the rotation, the presence of H, causes the mixing of the UV asymptotics in 

the metric functions to become slightly more involved. The rotation also causes the 

metric function expansions to become unwieldy at higher order, so we only perform 

the matching to 0(e~ 2p *) which is enough to render the ADM energy finite. Despite 

these technical complications the process proceeds as above. Namely, before taking the 

p* — > oo limit we fix coefficients so that the BPS and finite temperature metrics match 

at the boundary, we then take the p* — > oo limit and obtain a finite energy given by 



iaft 



5c 2 + C 2 
96tt 4 



(5.15) 



Thus, the energy remains the same as before the rotation. 

The above analysis was done for flavored finite-temperature backgrounds with 
s ^ 0. However, since the results do not depend on s they should also hold for 
backgrounds without flavor i.e for non-extremal deformations that, after the rotation, 
do have Klebanov-Strassler asymptotics. We have explicitly checked that this is indeed 
the case. 

In Fig. (19) we plot de/ds 14 and T as functions of the entropy density (or equiva- 
lently the "radius" of the horizon) for solutions before the rotation. The overlap of the 
two curves show that the numerical solutions satisfy the first law of thermodynamics 

de = Tds. (5.16) 



s denotes the entropy density S/V3 and should not be confused with the flavor parameter s. 
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Figure 19. deb e f/dsb e f (red line) and Th (blue dots) plotted vs. S& e /, for the flavored case 
with c+ = 50, s = 1. 



To verify that this is also the case after the rotation recall that we have already shown 
that the temperature and energy density remain the same (5. 14), (5. 15), (3. 6). It is easy 
to see that the entropy density is also unaffected by the rotation, 
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(5.17) 

5bef (5.18) 



p=ph 



where the last equality in equation (5.18) follows from the fact that e~ 8x vanishes at 
the horizon. 

Thus, the thermodynamical quantities e,s and T are not modified by the rotation 
and the first law will still hold. 

As discussed in section 3.4, our numerics associate a horizon temperature with 
a choice of c + and C<i- In Figs. (20) and (21) we plot the energy density versus the 
associated temperature for a solution with s — 1. The slope of these plots is the 
specific heat C v and is seen to be negative both before and after the rotation. Note 
that the fact that the first law is satisfied implies that we would have gotten the same 
answer had we chosen to use holographic renormalization methods. 
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The backreaction of the flavor branes does not play any role in the stability analysis 
as the temperature's dependence on s is essentially negligible compared to its depen- 
dence on C*2 and neither the energy nor the entropy depend on s. In Figs. (22) and (23) 
we show the energy density versus the entropy density S/V3. Note that in Figs. (22) 
and (23) the behavior of the energy vs. entropy is linear for large values of S (or equiv- 
alently for large values of the horizon "radius"). From the first law of thermodynamics, 
dE = TdS, we expect this to happen when the temperature is constant. 

Thus, we have shown that, despite of having very different asymptotics, the back- 
grounds after the rotation have thermodynamical properties similar to the D5 wrapped 
branes system. 




Figure 20. ADM energy density versus hori- Figure 21. ADM energy density versus hori- 
zon temperature (s = 1, c + = 50) zon temperature (s = 1, c+ = 3) 




Figure 22. ADM energy density versus en- Figure 23. ADM energy density versus en- 
tropy density (s = 1, c+ = 50) tropy density (s = 1, c + = 3) 
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6 Conclusions 



Let us summarize what we have done in this paper and point out possible future 
directions. 

• We found finite temperature solutions describing N c D5 branes wrapped on the 
S 2 of the resolved conifold. Unlike the generic wrapped D5 branes backgrounds 
the ones presented here have a dilaton that does not blow up at infinity but 
stabilizes to a finite value. 

• We identify a set of 11 UV and 6 IR parameters that determine the asymptotics 
of the solutions to any order. 

• It was not our goal to explore all the possible families of solutions in the 17 pa- 
rameter space but to restrict ourselves to those that have a UV behavior similar 
to the ones in [16]. Therefore, we fixed several of the UV parameters to appropri- 
ate values. We then solved the EOM's numerically using these UV expansions as 
boundary conditions and demanding regularity at the horizon. We also imposed 
that the constraint coming from reparametrization invariance is satisfied, to order 
1CT 6 or better, throughout all the interval. 

• Using U-dualities and these backgrounds as "seeds" we generate solutions with 
.D3 and D5 charge. To avoid issues related to the uplift of the smeared flavor 
branes, we show explicitly that after the U-dualities are applied the background 
obtained is a solution of the EOM's of D3 and D5 brane sources. 

• In the absence of temperature the backgrounds obtained are dual to interesting 
field theories that exhibit Seiberg dualities, Higgsing and confinement [16]. It is 
tempting to think of the non-extremal solutions found here as dual to finite tem- 
perature versions of the field theories in [16]. However, for non-zero temperature 
we have to first study the thermodynamical stability. We proceed to do that and 
find that the specific heat is negative and thus, they are unstable. 

• The solutions in [21] are a particular case (no flavor, a(p)=0) of ours, thus the 
thermodynamical instability found here applies also to [21]. 

• We show that our numerical solutions satisfy the First Law of Thermodynamics. 

• We find that the temperature grows as the non-extremality parameter decreases. 
Also, as we increase the area of the horizon the temperature goes to a small non- 
zero value. A similar behavior was found non-extremal MN [25] and indicates 



- 33 - 



the possibility of a phase transition. It is interesting to pursue this issue further 
and explore if the dualities can be used to map a phase transition in a wrapped 
fivebrane background to one in non-extremal KS [22]. Note that the tempera- 
ture dependence of Frg\ and the warp factor necessary to describe the chirally 
symmetric phase, T » T c , [27] emerges naturally in this framework due to the 
dualities. 

• It would be interesting to study wether there is a connection between our results 
and the ones existing in the literature regarding non-extremal fractional branes 
at an orbifold point [36]. 

• In [37] it was found that the lowest quasinormal mode in non-extremal Maldacena- 
Nunez is tachyonic. In this paper we have studied more general backgrounds than 
the one in [3 7]. The negative specific heat seems to indicate that a tachyonic mode 
is also present in these non-extremal flavored backgrounds, it would be interesting 
to confirm this. 

In [16], [17], [18], [19] it was shown that solutions of type IIB supergravity with 
sources generalize the Klebanov-Strassler baryonic branch; the dual field theory exhibits 
Seiberg dualities, Higgsing and confinement at different scales and was conjectured to 
describe a mesonic branch. One of the motivations of the present work was to find grav- 
ity duals of these interesting field theories. We find that starting from non-extremal 
wrapped D5 branes and applying the same dualities used to generate the extremal 
backgrounds produces thermodynamically unstable backgrounds. What is then the 
finite temperature gravity dual of the field theories found in [16], [17], [18], [19]? This 
fascinating question will undoubtedly lead to interesting physics and remains open. 
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A Appendix: The equations of motion 



In this appendix we briefly review the metric ansatz and solution method of [25], which 
we use to derive the Einstein equations and constraint given in section 2.2 for the fla- 
vored non-extremal backgrounds before rotation. 

The metric ansatz in [25], in Einstein frame, reads 

ds] IB = -Y x dt 2 + Y 2 (dx\ + dx\ + dx\) + Y 3 dp 2 + Y 4 (d6 2 + sin 2 6d(p 2 ) 
+ Y 5 (wi + adO) 2 + (u 2 -asinOdip) 2 ) + Y 6 (u 3 + cosOdip) 2 . 

The F( 3 ) gauge field ansatz is given by 
N, 



F (3) = t 



— (u\ + b d6) A (o>2 — b sin 8d(p) A (u 3 + cos 8d(f) 



N f 



+ b'dr A (-d6 A u x + sin Bd(p A Uq) + 11 — 6 — — r sin 6»d6» A^Au 3 



(A.1) 



where the Y^ a, and b are functions of p only. 

The complete action for type IIB supergravity with smeared flavor branes is then 

S 5grav 'S'sources; (-^■•2) 

with S'grav and Sources given by equations (2.10) and (2.12). 

We plug the metric and F( 3 ) ansatz into this action and integrate over all coor- 
dinates except p, drop the overall volume factor, and end up with a one-dimensional 
action with Lagrangian 

L = G ij( Y ) Y i Y j - U(Y) =T-U. (A.3) 
We can express the Y> in terms of nine other functions to make Gij diagonal: 

V — C 2z-Gx v — c 2z+2x V — c l0y-2z+2l 

i\ — e , r 2 — e , r 3 — e , 

F 4 = e 2^-2,+2 P +2 g) y 5 = e 2y-2z+2p-2 q ^ y & = e 2y-2z-8p^ ( A4 ) 

Y 7 = a, Y 8 = b, Y 9 = 4>. 
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We end up with the following expressions: 



T = e~ l (V 2 - 3x' 2 - 2z' 2 - hp' 2 -q 2 ~\ e~ Aq a' 2 - ^ e ^-^v h '2 _ 

Sy | e -12 P ^4q + g-4^2 _ ^2 + 2g 2^ _ e 10 P -2 g )2] _ ^-2p ^ 2g J 



16 1 V iV c 



^ + 2 (a-6) 2 | + Ar c e*/ 2 + 22+ ^- 4 P 



(A.5) 



Note that the N c dependence can be absorbed by shifting the dilaton, which we will 
do from now on. We will also write the ratio of the number of flavor branes to color 
branes as s = Nf/N c . Since I has no kinetic term, it is a pure gauge degree of freedom 
reflecting the remaining reparametrization invariance. Varying with respect to I one 
can set it to any value. 

We want to make connection with the metric ansatz given in equation (2.16) in Einstein 
frame, i.e. one in which 

Y x = e^ 2 e- 8x , Y 2 = e^ 2 , Y 3 = e ^ 2 e 2k e 8x , 
F 4 = e ^ 2 e 2h , Y 5 = (e^ 2 e 2 9)/4, Y 6 = {e^ 2 e 2k )/A, (A.6) 
Y 7 = a, Y 8 = 6, Y 9 = 0. 

In order to do this, we must pick the gauge I = — 4p — 4y + 4x + log 2, and also make 
the transformations 

y = }-(4 g + 4h + 2k - 10a; + 50 - 3 log 4), p = ^(g + h - 2k + log 2), 
q=-{-g + h + \og2), z = -^{-Ax + 4>) (A.7) 



After these transformations, the metric ansatz takes the desired form. In practice 
we first compute the equations of motion from the Lagrangian in the p, q, y, z, I, x 
variables, and then apply the above transformations to get equations of motion in 
the k,g,h,x,(p variables. This results in equations (2.21 - 2.28) for the EOMs and 
constraint, given in section 2.2. 

We note that the equations of motion and constraint are invariant under constant 
shifts of the dilaton 

+ (A.8) 

as well as shifts in the radial coordinate 

r — > r + tq. (A. 9) 
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B Appendix: UV asymptotics 
B.l Finite temperature UV asymptotics 

In the UV, we express the metric and gauge functions as a series in powers of e _2p / 3 : 
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i=l j=0 

The coefficients h it j, ■ ■ ■ , 6y are not all free, and inserting these expansions into 
the EOMs and constraint equation (2.21 - 2.28) allows us to express them in terms of 
11 independent coefficients, which we list here along with their order in the expansion: 



xi, (const), a^o (e" 8p/3 ), /i, (const), / 3 ,o (e" 8p/3 ), 

a 2 ,o (e~ 2p/3 ), a 4 , (e" 2 '), &4,o (e~ 2p ), 

ko,o (e 4p/3 ), k w (e- 8 ^ 3 ), 

h lfi (const), hx,\ (p). (B.2) 

Here we present expressions for the remaining non-zero dependent coefficients in 
the UV expansions, up to O (poly(p) x e~ 10p//3 ) . In order to simplify some of the higher 
order coefficients in what follows, we will assume that xxq = 1 and 02,0 = 0. 1 



.£>(e 4 "/ 3 ): 
»C(const): 



3 3 
^0,0 = 0^0,0 #0,0 = 7:^0,0 (B-3) 

o Z 



s 9s 

h,o = 2' 01,0 = -J - 4/l i,o, 9i,i = -4/ii,i, (B.4) 



15 The latter condition is required for a supersymmetric solution, and we set it equal to zero also in 
our finite temperature solutions, such that C2 is the only parameter deforming the BPS backgrounds. 
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B.2 Supersymmetric UV asymptotics 

Here we summarize the BPS, zero-temperature UV asymptotics (up to 0(e~ Ap / 3 )) of the 
backgrounds corresponding to D5s wrapped on the S 2 of the resolved conifold, with the 
addition of smeared D5 flavor sources, as described in [16, 28]. These serve as matching 
conditions for our finite temperature UV asymptotics. In addition, they furnish us with 
a family of reference backgrounds that are used in the energy calculations of section 5. 

The BPS asymptotics are given in terms of the free variables Q ,c + ,c^, p and 
/i,o after performing a UV expansion of e.g. equations (3.5) of [16]. 16 The finite 
temperature asymptotics of section 3.1 are a deformation of these BPS asymptotics 



16 In [16], the parameter <p is used instead of /i,o- We remove the c+ dependence of the asymptoti- 
cally constant part of the dilaton by setting e A ^° = 2/i,qc5_/3. See equation (B.46) of [16]. 



- 40 - 



with p = 0, Qq = —N c + Nj/2 and /i,o = 1- In the energy calculations of section 5, 
we also set po = 0, but Q , c+, c_ and / i>0 are left free to allow for the matching. 

What follows, then, is the form of the UV BPS asymptotics with p = 0, but with 
Qo, c + , c_ and /i 5 o left free. 17 
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(B.9) 



C Appendix: Comments on numerical procedure 

The UV-shot solutions of section 3.3 were obtained by using Mathematica's NDSolve, 
with UV boundary conditions determined by the expansions given in equation (B.l) up 
to 0(e- 8p ) for e 8x , C( e - 26 ^ 3 ) for a and 6, C( e - 20 ^ 3 ) for e 2k ,e 29 and e 2 \ and 0(e' 8p ) 
for e 4<?i . Using an expansion taken out to this order, NDSolve at WorkingPrecision = 
70 finds a horizon radius ph which is fairly independent of the value of p^ chosen (~ 
1 part in 10 -5 as p^ is varied from 6 to 9). The constraint equation (2.28) is of order 
1CT 8 when evaluated on a typical UV-shot solution, but climbs near the horizon; on 
the horizon shot solution the constraint is of order 10 -14 with a maximum of ~ 1CT 6 
at the horizon — see Figs. (24) and (25). The degree to which the constraint equation 
is violated near the horizon also depends on the size of Ci. This stems from the fact 
that many higher order terms in the UV expansion are proportional to C*2, so any finite 
order truncation of the UV asymptotics becomes less accurate for larger C2 values. 

17 c_ appears only at 0(e -8p / 3 ) and higher in the metric functions, which we omit here for brevity. 
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When matching our UV-shot solutions to a series expansion near the horizon, we 
choose to exclude the a and b functions, the dilaton, and their derivatives from the 
mismatch function. The reason for excluding the a and b functions is that unlike the 
metric and dilaton functions — whose UV and horizon shot solutions agree quite well 
away from the horizon in what Minimize considers the best matched case — a and 6's 
UV and horizon shot solutions differ over the entire interval in the best matched case, 
as shown in Fig. (26) (although a and b are typically small in magnitude compared 
to x,g,h,k and 0). The UV-shot dilaton also displays divergent behavior near the 
horizon, and drastically decreases the match's accuracy if included. Including it in the 
match also produces a horizon-shot dilaton which stabilizes to a UV value different than 
e Ar ^\ Pac ~ 1. We remedy this by using the invariance of the EOMs under a constant shift 
in the dilaton. After the matching is performed, we pick an that gives a horizon-shot 
dilaton that agrees with the UV value of e 4 ^| Poo ~ 1. For the solutions we have found, 
including a, b and the dilaton allows for the mismatch to be minimized to only ~ 10~ 2 , 
and the resulting horizon shot solutions for (x, g,h,k, f) differ significantly from the 
UV shot solutions. Excluding a, b and <j) from the mismatch function, Minimize 
is able to reduce the mismatch to ~ 1CT 7 , and the resulting horizon shot solutions 
for (x, g, h, k, <p, a) agree very closely with the UV shot solutions. This will still leave 
us with a horizon-shot b that tends to exponentially diverging behavior at large p. 
However, b can be tuned to eliminate this divergence, with negligible effect on the 
behavior of the other functions (x,g, h, k, f) — see Fig. (26). 

Performing our matching close to the horizon (e.g. e = .15) necessarily involves 
greater errors in the UV-shot solutions. This shows up as horizon-shot solutions that 
disagree with the UV-shot solutions near p^. To remedy this, we instead begin by 
matching at a higher e = .7 — this eliminates the disagreement near p^. We then 
use the resulting values of the matched horizon parameters as seeds for a new match 
performed at a lower e = .15, where in addition we constrain the allowed variance of 
(xx, go, ho, k , f ) around the seed values until the resulting e = .15 match gives horizon 
shot solutions which agree with the UV shot solutions at p^. After doing this, we find 
that the high-e-seeded e = .15 match produces a mismatch of the same order (10 -7 ) as 
a non-seeded e = .15 match, with the added benefit that the high-e-seeded horizon-shot 
solutions are in good agreement with the UV-shot solutions at p^ 

The numerical method we use could undoubtedly be improved. However, as pre- 
sented here it is enough to identify the behavior of the horizon temperature as a function 
of the non-extremality parameter C<i, which is necessary to study the stability of the 
background. Other methods have been used, such as [24]. 
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Figure 24. Constraint evaluated on a horizon Figure 25. Constraint evaluated on same so- 
shot solution with c + = 50, C2 = 5000, s = 1. lution, near horizon at ph ~ 3.194024. 




Figure 26. Horizon-shot and UV-shot a and 6 after tuning bo. brown - UV-shot a, black - 
UV-shot b, blue - horizon shot a, red - horizon-shot b. s = 1, c = 3, C2 = 800000. 



D Appendix: U-Duality for the flavored, wrapped D5 branes 
black hole 

Here we describe in more detail the steps of the U-duality procedure as applied to 
backgrounds of the type mentioned in equations (2.16) and (2.15). We will consider 
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the rotation in string frame. We begin with 



dsi = e 



+ 



Hp) 

e 2k(p) 



- e- 8x(p) dt 2 + dxj + dx\ + dx\ 



+ ds 6 , 
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N f 



(D.l) 



b'dp A (-d9 Awi + sin 9dip A w 2 ) + (1 - b(p) 2 - —£) sin9d9 Ad<pAu 3 



(D.2) 



We will supress obvious p dependence to avoid cluttering. We begin by T-dualizing in 
the x±,X2,X3 directions, which results in the type IIA background 



d- s iiA — e 



- e- 8x dt 2 ^ + e-^{dx\ + dx\ + dxj) + ds 2 6 , 
F (6 ) = F (3 ) A dxi A dx 2 A dx 3 ->■ F (4) = e 2< V 4:r * 6 F (3) A dt. (D.3) 



e 2 ^ = e 
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(D.4) 



Now, we lift this to M-theory: 

ds 2 u = e- 2 * /3 dx 2 u + e<^ 3 [ - e-^e^dt 2 + e^(dx 2 + dx\ + dx 2 ) + ds 
G (4 ) = e~ 4 V* * fi A eft. 



6 -^(3) 



(D.5) 
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Boosting in the (t, xn) directions according to 

dt cosh f3dt — sinh (3dxu, dxn — > — sinh (3dt + cosh (3dxn, (D.6) 
we rewrite the boosted metric as 



ds^ 



,0/3 



e ^{dxl + dx\ + dx\) + dsl + Adt 2 + Bdi^ + Cdtdxi 



G 



(i) 



e- 4 V» * 6 F (3) 



cosh /3<it — sinh (3dxn 



(D.7) 



where 



A = e 4 * /3 [sinh 2 /3e" 2 * - e~ 8x cosh 2 = e 4 * /3 [cosh 2 /3e" 2 ^ - e~ 8x sinh 2 /3], 

C = -2cosh/3sinh/3e 4 * /3 [e~ 20 - e~ 8x }. (D.8) 

Before reducing to IIA, it is useful to write equation (D.7) as 

g tt dt 2 + B^e^ie^idxj + dx\ + dx\) + dsl)] + B(dx u + a t dt) 2 , 



dsl, = B- 1/2 



G (4) =e- 4 V* * 6 F (3) 
where we have defined 

a t = 



(cosh + a t sinh f3)dt — sinh (3(dxn + a t dt) 
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2B' — 4VB 
Now we reduce to IIA, obtaining in string frame, 
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(D.10) 
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^(2) — a [dp A (it. 

Next, we T-dualize back along the Xi,x 2 ,x 3 directions, and obtain 
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Finally we use the definitions for A, B, C and a t , and take the limit j3 — > oo. This is 
the field theory limit, where the warp factors vanish at infinity. We then rescale 



N = N cosh/3, Xi ->• ^cosh /3 V N a 1 Xi. (D.13) 

With all of the above, the — > oo limits are finite and the final solution is given by 
equations (4.1) and (4.2) after transforming to Einstein frame 18 . 

E Appendix: General form of 

In the extremal case [16], the SU(3) structure fixes the form of -8(2)- in our solutions 
5(2) should reduce to the one in [16] when T = 0. The following ansatz is compatible 
with that requirement: 

B{2) = bi(p)uJ 3 A dp + b 2 (p)e 1 A e 2 + h(p)e 1 Au 2 + h(p)e 2 Aui + h(p)u 1 A oo 2 . (E.l) 

We find it convenient to introduce a new function b 2 (p) and parametrize b 2 as, 

b 2 (p) = b 2 (p) + (1 - a{pf)b 5 {p) (E.2) 

Therefore, 

B{2) = &i(p)cJ 3 Adp+(6 2 (p) + (l-a(p) 2 )&5(p))eiAe2+63(p)eiAa;2+&4(p)e2Aa;i+65(p)wiAtI;2 

(E.3) 

With this notation it is easy to show that 

B i2) = B {2) \ b5=Q -d[h(p)iJ 3 ] (E.4) 

Thus, 65 is really a gauge choice that does not affect the value of H( 3 y 
From (4.1) we have, 



#(3) = -e- 4 ^ * 6 / {3) 



Demanding that if( 3 ) = d[B^} results in five equations, 
cat(0)(6 3 -64) = 
h + ab 5 - \e 2<t> ~* x b' = 



(E.5) 



2e 2<t>-2g+2h _ ^ + ^ + ^ 3 + h + 2ab5 y _ ^ _ _ a 2) = o 

e 20 (-a + b) - 2(a&i + & 5 a' + 63) = 

1 2*+2g-2h (_ 1 + s _ a 2 + 2fl6 ) + bi _ & / = o (E 6 ) 



18 



To avoid cluttering the notation, in equations (4.1) and (4.2) the rescaled quantity is called N c 
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Three of them are easily solved by setting, 
h = h 

h = - 1 - {8ab, - e 2 ^b>) 

o 

6 X = I e - 2/l+29+2< ^ (1 - s + a 2 - 2ab) + b' 5 = (E.7) 

we are then left with two equations, one of them is just the equation of motion for b(p). 
The other one is a differential equation for 62 (p), 

- e* x + 29+2h h{p) = (E.8) 

Note that equations (E.7-E.8) determine B2 in terms of an arbitrary function 65 that 
we have shown is just a gauge choice. In order to make contact with [16] we choose, 

b' 5 = -e 2 ^ 2h [Ae 2h+2k + e 29 (-l + s - a 2 + 2ab)) (E.9) 



which gives b\ 



\J2<t>+2k 
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